WHAT IS A SEQUENCE OF NILSSON TYPE? 



STAVROS GAROUFALIDIS 

Abstract. Sequences of Nilsson type appear in abundance in Algebraic Geometry, Enumerative Combina- 
torics, Mathematical Physics and Quantum Topology. We give an elementary introduction on this subject, 
including the definition of sequences of Nilsson type and the uniqueness, existence, and effective computation 
of their asymptotic expansion. 
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1. Sequences of Nilsson type: definition 

Sequences of Nilsson type are the ones that are asymptotic to power scries in powers of 1/n and logn. 
They appear in abundance Analysis (in asymptotic expansions of integrals), in Mathematical Physics and 
in Algebraic Geometry (in relation to the Gauss-Manin connection) ; see for example [An2[ IMli IM2[ IPhl ISa) . 
They also appear in Enumerative Combinatorics (see |FS[ |WZ2 ( Ga2j) and in Quantum Topology. For 
instance, the Witten-Reshetikhin-Turaev invariant of a closed 3-manifold is a sequence of complex numbers 
that depends on the level, and it is expected to be of Nilsson type; see |Wi[ IFG1 IGa31 IRol ILR1 1 AH) . In 
addition, the Kashaev invariant of a knot is expected to be a sequence of Nilsson type; see |KT[ IAH1 ICGj . 
The quantum spin network evaluation at a fixed root of unity is known to be a sequence of Nilsson type; 
see [GV1| IGV2] . For a general discussion of perturbative and non-perturbative invariants of knotted objects 
that are expected to be sequences of Nilsson type, see [Gal] . 

There is a close connection between sequences of Nilsson type and multivalued analytic functions with 
quasi- unipotent monodromy; see for example Theorem Q] below. 

Several people familiar with the ideas of Quantum Topology have asked for a self-contained definition of 
sequences of Nilsson type and their asymptotics, its uniqueness, existence and effective computation. 
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Asymptotics of sequences is a well-studied subject of analysis that goes back at least to Poincare; see 
for example [O] [C] IM1) . Since we could not find a reference for sequences of Nilsson type in the existing 
literature, we decided to write this introductory article. It concerns the asymptotic expansion of sequences 
which are relevant in Quantum Topology, and may serve as an elementary introduction to asymptotics. We 
claim no original results in this survey paper. 

In order to define sequences of Nilsson type, we need to introduce Nilsson monomials h u (n) indexed by a 
well-ordered set fi, and a finite set A of complex numbers of equal magnitude. 

For a natural number d € N, a finite subset S of the rational numbers consider the well-ordered set 
= (S + N) X {0, 1, . . . , d} indexed by (a, ff) < (a', f}') if and only if a < a' or a = a' and /3' < f3. SI has 
the order type of the natural numbers. In particular, for every u) G O, the set of elements strictly smaller 
than u> is finite. Consider the Sl-indexed family of monomials of Nilsson type given by: 

(1) K(n) = na 

for u) = (a,/3) £ S7. It is easy to see that \im n h u i(n)/h u (n) = (abbreviated by h u i (n) / 'h u (n) = o(l), and 
also by h u (n) 3> h U '(n)) if and only if u) < lo' . This and all limits below are taken when n goes to infinity. 

Fix a finite set A of nonzero complex numbers of magnitude r > 0. Let c Ut x be a collection of complex 
numbers indexed by SI x A. 

Definition 1.1. (a) With the above notation, for a complex-valued sequence (a n ) the expression 

(2) a„~ EM n )X>oa A " 

ujGSI AeA 

means that 

• for every uj G Q we have: 

(3) \a n r- n - £ Mn)X;<V.A(Ar- 1 ) n ] j-\t = o(l) 

\ u'<u AeA J u ^ > 

• Cui,\ 7^ for some (u>, A) € fl X A. 

(b) We say that a sequence (a n ) is of Nilsson type if there exist SI, A and Cu,\ such that ((2|) holds. 

We will say that an asymptotic expansion ([2]) is SI x A-minimal if 

• For every AeA there exists uj £ SI such that c w ,a 0. 

• For every u) € SI there exists A £ A such that c w ,a 0. 

By considering a subset of A or SI if necessary, it is easy to see that every asymptotic expansion has a minimal 
representative. 

2. Uniqueness 

Our first task is to show that a sequence of Nilsson type uniquely determines SI, A and the coefficients 
c u _\. The key idea is the following elementary lemma. 

Lemma 2.1. If A is a finite subset of the unit circle and 

(4) E caA" = o(l) 

AeA 

holds for some complex numbers c\, then c\ = for all A G A. 

Proof. Divide © by A™ for some Ai € A. Then we have ca x + X^a^Ai c a(A/Ai)" = o(l), where A/Ai 7^ 1. 
So, 



(<*!+£ C A(A/Ai) fe j =o(l) 

k=l \ A#A X / 
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1 ^ 1- A/A! " +i 
c *i + nS> /-A/A, =°W- 



< C 



By averaging, it follows that 

1 - (A/Ai)" +1 
A/A 

Thus, = 0. Since Ai was an arbitrary element of A, the result follows. □ 
Lemma 2.2. // (a n ) satisfies @ then 

(5) limsupKI 1 /" = r. 

n 

Proof. Since c Wj a 7^ for some (to, A) e x A, without loss of generality assume that c WOi a 7^ for some 
AeA where ujq is the smallest element of tt. Equation ([2]) for lj = ujq gives that 

(6) ^_^^ iA (Ar-^ = o(l) 

ASA 

Now Ar _1 are on the unit circle. It follows that 

huo (n) 

for some C > 0. Since lim n /^(ri) 1 /™ = 1 for all u> € O, it follows that 

limsup lonl 1 /™ < r 

n 

Since some c^ 0! a is nonzero and Ar _1 are on the unit circle, Lemma[2J]implies that lim„ J2\eA c w ,A(Ar -1 )" 7^ 
0. Since the sequence is bounded, it follows that there exists a subsequence n& such that 

limV^jtAr-'r^'/O 

AeA 

Combined with Equation (j6]), it follows that 

it 

The result follows. □ 

In particular, Lemma POl implies that sequences of Nilsson type satisfy limsup n la,^ 1 /™ > 0. 
Proposition 2.3. Suppose that 

(7) a n ~ ^Wj]^" 

wen AeA 

and 

(8) E m») E c -',a' a '" 

u'eff A'eA' 

are f2 x A-minimal and £1' x A' -minimal asymptotic expansions. Then Q = Q! , A = A'. Moreover, for all 
(ui, A) € fl x A we have c^^x = d u x . 

Proof. Let wo and lu' denote the smallest elements of ft and f2'. Lemma [2?2l implies that r — r' where r and 
r' are the magnitudes of the elements of A and A' respectively. Equation ([3]) for luq and uj' implies that 

(9) ^ ~ E ^r-r = o(l), ^1 _ ^ 4, >A ,(AV-)« = 
ft(JoW AeA < [ > A'eA' 

If wo 7^ we may assume that ujq < uj' q . In that case, observe that h u ' (n)/h uo (n) = o(l). Multiply the 
second equation above by h u > (n)/h Uo (n) and subtract from the first. It follows that 

-E c -A(Ar-T + ^ E c' u ^'r-r = o(l) 
AeA u ° y ' A'eA' 
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Since h u ' o (n)/h Uo (n) = o(l), it follows that 

X)^.A(Ar- 1 )" = o(l) 
AeA 

Lemma 12.11 implies that c Wo ,a = for all A contrary to our minimality assumption of (jTJ) . It follows that 
uq = uj'q. Subtracting, Equation (|9]) implies that 

- y Zc UOlX {\r- i r+ J2 ^ M x ' r ~ 1 ) n = °( 1 ) 

ASA A'eA' 

Lemma |2. II implies that if c UOt \ ^ for some A € A, then A € A' and moreover c UOt x = c' LJg x , . 

An easy induction on uj <E f2 proves the following statement. For every w G O, the following holds. If 
c u ,A 7^ for some A £ A, then AeA' and ui € ft' and c w ,a = A . 

The minimality assumption and the above statement implies that Q — tt' and A = A' and cx, u = c' x for 
all (uj, X) £ ft x A. □ 



Remark 2.4. Proposition 12.31 proves uniqueness in a non-effective way. We will come back to the problem of 
computing c Wj a later on. 

3. Alternative expression for sequences of Nilsson type 
If (a n ) is a sequence of Nilsson type, we can write ^ in the form: 

(10) a„~ Y, ^(lognfSx^gx^A 1 /™) 

X,a,0 

where 

• the summation in (1101) is over a finite set, 



• the growth rates A are complex numbers numbers of equal magnitude, 

• the exponents a are rational numbers and the nilpotency exponents j3 are natural numbers, 

• the Stokes constants S\ >a ,p are complex numbers, 

• Q\,a,p{%) G C[[a;]] are formal power series in x with complex coefficients and leading term 1. 

Remark 3.1. In the definition of a sequence of Nilsson type, we may additionally require that 

• A is a set of algebraic numbers, 

• the formal power series gx,a,p{x) is Gevrey-1, i.e., that the coefficient of x k in gx, a , p(x) is bounded 
by C k kl for all k, where C depends on gx,a,/3, 

• the coefficients of the formal power series gx,a,p(x) lie in a fixed number field K , 

These additional requirements hold for the evaluations of classical spin networks, see [GVlj . as well as 
Sections 0] and 16.11 below. 

Example 3.2. For example, if d = 1 and S = {1,3/2}, then Q = (1 + N) U (3/2 + N) and we have: 

log n 1 log n 1 log n 1 

— — > - > > — — > — %r- > — >... 

If A = {k, fi, v\, the asymptotic expansion (|T0|) of a sequence of Nilsson type becomes: 

logn^ /1\ log (I 

AeA x ' AeA x 



AeA x ' AeA x 

where g\.j(x) € C[[x]] are formal power series in x and Sx,j are complex numbers 
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4. Existence 



In this section we will prove that a sequence is of Nilsson type, under some analytic continuation assump- 
tions of its generating series. This is a well-known argument (see for example, [Cl ICG1 iFSl IGM1 IGIKM1 IMlj ) 
that consists of the following parts: 

• apply Cauchy's theorem to give an integral representation of the sequence, 

• deform the contour of integration to localize the computation near the singularities of the generating 
function, 

• analyse the local computation using the local monodromy assumption of the generating function. 

Let us give the details of the above existence proof. Since sequences of Nilsson type are exponentially bounded 
(as follows from Lemma l2.2p . fix an exponentially bounded sequence (a n ) and consider its generating series 



G{z) is an analytic function for all complex numbers z that satisfy \z\ < 1/R. Suppose now that G has 
analytic continuation on a disk of radius r with singularities at finitely many points re, A, fj,, u, . . . . Suppose 
also that G has further analytic continuation on a disk of radius r+e minus finitely many segments emanating 
from the singularities radially as in the following figure. 



Assume in addition that G has quasi-unipotent local monodromy at each singularity A, \i, v, k on the circle of 
radius r (i.e., the eigenvalues of the local monodromy are complex roots of unity). 

Theorem 1. Under the above assumptions, the sequence (a n ) is of Nilsson type. 

Corollary 4.1. Suppose that G{z) = X^^Lo a « z ™ * s a multivalued analytic function on C \ A (where A C C 
is a finite set) which is regular at z = 0, and has quasi-unipotent local monodromy. Then, (a„) is a sequence 
of Nilsson type. 

Remark 4.2. We know of at least three ways to show that a germ G(z) of an analytic function can be 
analytically continued to the complex plane, namely 

(a) G satisfies a linear differential equation, see for example |Ga2i Thm.l] reviewed in Section HTT1 below. 
For examples that come from Quantum Topology (specifically, spin networks) see |GV11 IGV2] . 

(b) G satisfies a nonlinear differential equation. See for example the instanton solutions of Painleve I 
studied in detail in IGIKM] and the matrix models of [GMj . 

(c) G is resurgent. See for example the Kontsevich-Zagier series studied in detail in |CGj . and more gen- 
erally the arithmetic resurgence conjecture of [Gal] for sequences that appear in Quantum Topology. 

Proof, (of Theorem[T]) We begin by applying Cauchy's theorem to give an integral representation of (<z n ). If 
7 is a circle of radius less than 1 / R that contains the origin, then we have: 



oo 



(11) 




71=0 




(12) 
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We can deform 7 to a contour C which consists of a contour T~i\ around each singularity A and finitely many 
arcs 7r+ e of the circle of radius l/(r + e) as in the following figure. 



(13) 





\ 




c /"""l 








y 
















\ ^ M 
\ r 




7tiu J 









The contours 7-L\ are known as Hankel contours in Analysis (see [O]) and Lefschetz thimbles in Algebraic 
Geometry (see [Phi ISa] ). Cauchy's theorem implies that 



(14) 



— V 

2ni ^ 

xeA 



G{z) , 1 



7r+ e 



G(£) 
7 n+l 



dz 



The above expression is exact, and decomposes the sequence (a n ) into a finite sum of sequences (one per 
singularity) and an extra term. Of course, there is nothing canonical about this decomposition, since the 
size of the Hankel contour and 7 r + £ depends on e. One could make the decomposition nearly canonical by 
using Hankel contours that extend to infinity, but even so there are choices of directions to be made, and we 
will not use them here. 

The integral over j r ^ e can be estimated by 0((r + e)~ n ) since G is uniformly bounded on the arcs 7 r . £ . 
Since we assume that the local monodromy of G(z) around a singularity is quasi-unipotent, it follows (see 
[Ml] ) that modulo germs of holomorphic functions at zero, G(A + z) has a local expansion of the form 

G(\ + z)= z a \\ogzfh a .^{z) 

a',p> 

where the summation is over a finite set, a' G Q, f3' G N and h a >^>(z) are germs of functions analytic at 
z = 0. For a germ f(z) of a multi-valued analytic function at z = 0, let Ao/ denote its variation defined by 

A f(z) = lim /(e 2 " e z) - lim f(e 2 ™z) 

e-vl £— >0 

(see [Ml] ) when z is restricted on a line segment [0, e). The variation of the building blocks z a and (log 2)^ 
are given by 

!(e 27rm -l)z Q aeQ\Z 
S aeZ\N, Ao(logz) = 2m 

a G N 

where So is the Dirac delta function (really, a distribution). For a singularity A of G(z), let A\G(z) denote 
the variation of G(X + z). It follows that for z in the line segment of Figure (p~3|) . we have 



(15) 



A x G(z) 



a,0 



z Q (logz) f 



fc-l 



k=0 



where the sum is over a finite set {a, /?}, qgQ\M, /ieN and X^fclo Sfc^o c q a k zk are g erms 01 analytic 
functions at z = 0. When a G Z\N, we can deform the Hankel contour into a small circle centered around A 



WHAT IS A SEQUENCE OF NILSSON TYPE? 



7 



and apply Cauchy's theorem. For the remaining cases a € Q \ Z, a change of variables z A(l + z) centers 
the Hankel contour at zero and implies that 



(16) / -^-jdz = A~" / v " dz = \- n / r . \ ' dz 

A Beta-integral calculation gives that 



o (l + z) n+1 Jo (1 + *) 7 

r( 7 )r(n + i- 7 ) 



-riz 



(l + z)"+! r(n + l) 



and therefore 



f - rwrt^ i - 7) (1 + 0((f + 

Jo (l + z)™+ 1 r(n + l) 

More generally, for a natural number /3 let us define 



/>oo 7 — 1 

Then, we have 

(18) ^(n) = r(7) ^^- 7 W ») 

where pp(^f,n) is a polynomial in the variables ip( k \n + 1 — 7) and ^^(7) with rational coefficients, where 
^(2) = r'(z)/r(z) is the logarithmic derivative of the r-function. For example, we have: 

p (n) = 1 

Pi(n) = --0(n+l- 7 ) + ^(7) 

p 2 (n) = V(« + l-7) 2 +V j(1) (« + 1 -7)-2'0(7)'0(« + l-7) + ^(7) 2 + ^ (1) (7) 
Compare also with [Mil Eqn.I.4.2] and M2, Eqn.7.5]. What is important is not the exact evaluation of 
If, pin) given in (TT51) . but the fact that the sequence 7 7) ^(n) is of Nilsson type. This follows from the fact 
that we have an asymptotic expansion (see [O]): 

£(n + 1-7) 1 A 7 2 - 7 3 7 4 - 2 7 3 - 3 7 2 + 2 7 
1 j r(n + 1) ~ nT V 2n 24n 2 

Alternatively, one may show that the sequence 7 7i( g(n) is of Nilsson type by a change of variables z = e t — 1 
which gives 



where ^ ^ 

A 7 ,^(t) = (^-^) (log - logt ) di 

is a function which can be expanded into a polynomial of log t with coefficients functions which are analytic 
at t — 0. Expand ^4 7ji g(t) into power series at t = and interchange summation and integration by applying 
Watson's lemma (see [O]) to conclude that J 7)/ g(n) is of Nilsson type. 

Replace AxG(Xz) by (|15p in ([T6")) and interchange summation and integration by applying Watson's lemma 
(see [Q]). It follows that 



55 £ ^~*-"E>s»> s X>«4 

J Ha a,/9 fc=0 

Equation f| 14f) cocnludes that (a n ) is of Nilsson type. Strictly speaking, the above analysis works only when 
$t(a) > — 1. This is a local integrability assumption of the Beta-integral. The asymptotic expansion @ 
remains valid as stated even when 3?(a) < — 1 as follows by first integrating G(z) a sufficient number of times, 
and then applying the analysis. This is exactly what was done in |CG1 Sec. 7] at the cost of complicating the 
notation. □ 
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Remark 4.3. Since the sequence (c' a p A k ) as a function of k is exponentially bounded and the asymptotic 
expansion (IT^l) is Gevrey-1, it follows that the sequence {ca,p,x,kj is Gevrey-1. Moreover, if the sequence 
[d a a A k ) lies in a number field K, then we can write the asymptotic expansion of («„) in the form (fT0|) 
where S a p t x are polynomials (with rational coefficients) of values of logarithmic derivatives of the Gamma 
function at rational numbers. 



5. G-FUNCTIONS 

5.1. G-functions: examples of sequences of Nilsson type. In |Ga2. Thm.l] it was proven that balanced 
multisum sequences (which appear in abundance in Enumerative Combinatorics) are sequences of Nilsson 
type. The proof uses the theory of G-functions which verifies that the generating scries of balanced multisum 
sequences satisfies the hypothesis of Corollarv l4.ll Let us give the definition of a balanced multisum sequence, 
a G- function and an example. 

Definition 5.1. (a) A term t n .k in variables (n, k) where k = [k\, . . . , k r ) is an expression of the form: 

r ,7 

(20) V^CjJJCf* JjA^n,*;)^ 

»=i j=i 

where C; G Q for i = 0, . . . , r, ej — ±1 for j = 1, . . . , J, and Aj are integral linear forms in the variables 
(n,k) such that for every n € N, the set 

(21) supp(t„,.) := {keZ r \ Aj(n, k) > 0, j = 1, . . . , J} 

is finite. We will call a term balanced if in addition it satisfies the balance condition: 

J 

(22) E e ^=°- 

.7=1 

(b) A (balanced) multisum sequence (a n ) is a sequence of complex numbers of the form 

(23) a n = W 

fcgsupp(t„,.) 

where t is a (balanced) term and the sum is over a finite set that depends on t. 

For example, the following sequence is a balanced multisum 

/ n \ 2 {n + k\ 3 {n + l\ _ v (n + k)l 3 (n + l)\ 

[ ' " ^\k + l) \ k ) \ I ) k\Hln\ 2 (k + iy. 2 (n~ k - 

where the summation is over the set of pairs of integers (k, I) that satisfy < k, I and k + I < n. 

Let us now recall what us a G-function. The latter were introduced by Siegel in [Si] with motivation 
being arithmetic problems in elliptic integrals, and transcendence problems in number theory. For further 
information about G-functions and their properties, see |Anl| [A"n2] . 

Definition 5.2. We say that series G(z) = J^^Lo a nZ n is a G-function if 

(a) the coefficients a n are algebraic numbers, 

(b) there exists a constant G > so that for every n £ N the absolute value of every conjugate of a n is 
less than or equal to G", 

(c) the common denominator of oq, . . . ,a n is less than or equal to G™, 

(d) G(z) is holonomic, i.e., it satisfies a linear differential equation with coefficients polynomials in z. 

G-functions satisfy the hypothesis of Corollarv l4.ll see |Anl[|A"n2] . Indeed, they satisfy a linear differential 
equation which analytically continues them in the complex plane. Moreover, the arithmetic hypothesis 
ensures that the local monodromy is quasi-unipotent. We can now state the main result of |Ga2j . 
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Theorem 2. [Ga2] (a) If (a n ) is a balanced multisum sequence, its generating function G(z) = X^^Lo a « z ™ 
is a G-function. 

(b) In that case, it follows that (a n ) is a sequence of Nilsson type. 

The reader may have noticed that we defined the notion of a sequence of Nilsson type only when 
limsup Icinl 1 /' 1 > 0. In case the generating series is a G-function, the remaining case is taken care by 
the following lemma. 

Lemma 5.3. If G(z) = X^^Lo a ™ z ' 1 * s a G-function and limsup la,^ 1 /" = 0, then a n — for all but finitely 
many n. 

Proof. The assumption implies that G(z) is an entire G-function. Since those are regular-singular at infinity, 
it follows that G(z) is a polynomial; see also |Anll lA"n"2] . The result follows. □ 

5.2. Classical spin networks: examples of G- functions. In |GVl[[GV2j it was proven that the evalua- 
tion of a quantum spin network at a fixed root of unity is a balanced multisum sequence, and consequently 
it is a sequence of Nilsson type. 

6. Effective computations 

6.1. Exact computations. Proposition ^. 31 is a uniqueness statement about the asymptotics of a sequence 
of Nilsson type, and Theorem Q] is an existence statement which is not effective. There are two types of 
effective computations, exact and numerical. The exact computations use as an input a linear recursion 
relation of the sequence. The following proposition is elementary and is discussed in detail for example in 
[FS1IWZ2] . 

Proposition 6.1. Given a linear recursion relation for a sequence (a n ) of Nilsson type, one can compute 
exactly \,a,0 and the power series g a ,b.\{x) that appear in Equation (I10|) . 

In particular, a linear recursion relation computes exactly the asymptotics of a sequence of Nilsson type, 
up to a finite number of unknown Stokes constants. 

To apply Proposition 16.11 one needs to find a linear recursion for a sequence (a„). This comes from the 
fundamental theorem of Wilf-Zeilberger which states that a balanced multisum sequence is holonomic, i.e., 
satisfies a linear recursion with coefficients polynomials in n; see [Zj IWZ11 IPWZj . The proof of the above 
theorem has been computer implemented and works well for single sums and reasonably well for double 
sums; see [PWZ1 iPRll IPR2] . As an example, consider the following sequence from [GVll Sec. 10] 

n! 6 ^ (-l) fc (fc + l)! 
a ™ ~ (3n + l)! 2 ^ {k - 3n)! 4 (4n - fc)! 3 

Using the language of |GVlj . (a n ) is the evaluation of the tetrahedron spin network (also known as 6j-symbol) 
when all edges are equal to n. The command 

<< zb.m 

loads the package of |PR2] into Mathematica. The command 

teucl[n_, k_] :=n! A 6 / (3n + l) ! A 2 (-l) A k (k + 1) ! 1/ ((4n-k) ! A 3 (k-3n) ! A 4) 

defines the summand of the sequence (a„), and the command 

Zb[teucl[n, k] , {k, 3 n, 4 n} , n, 2] 

computes the following second order linear recursion relation for the sequence (a n ) 

-9 (1 + n) (2 + 3n) 2 (4 + 3n) 2 (451 + 460n + 115n 2 ) a[n] + 

(3 + 2n) (319212 + 1427658n + 2578232 n 2 + 2423109 n 3 + 1255139 n 4 + 340515 n 5 + 37835 n 6 ) a[l+n] - 
9 (2+n) (5 + 3n) 2 (7 + 3n) 2 (l06 + 230n + 115n 2 ) a[2 + n] ==0 

This linear recursion has two formal power series solutions of the form 
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1_ ( -432±31i%/2 109847 T 22320iV2 -18649008 ± 4914305iv / 2 
a± ' n ~ ± \ + 576^ + 331776n 2 + 573308928n 3 

14721750481 ± 45578388960^ -83614134803760 ± 7532932167923^ 

660451885056n 4 + 380420285792256n 5 

-31784729861796581 =p 212040612888146640^ 



where 



657366253849018368n 6 



329T460V2 _ Tl6arccos(1/3) 
l± 729 




are two complex numbers of absolute value 1. The coefficients of the formal power series a± :U are in the 
number field K = Q{V^)- 

6.2. Numerical computations. When a sequence (a„) is given by a multi-dimensional balanced sum, the 
computed implemented WZ method may not terminate. In that case, one may develop numerical methods 
for finding A, a, f3 as in Equation (jlOp . An example of this method is the asymptotics of the evaluation of the 
Cube Spin Network that appears in the Appendix of |GV1| . Effective methods for numerical computations 
of asymptotics have been developed by several authors, and have also been studied by Zagier. 
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